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1.  Introduction 

Although  high-power  fiber  amplifiers  are  commereially  available  at  the  lOkW  level,  their 
use  in  direeted-energy  weapon  (DEW)  systems  require  narrow  linewidth,  whieh  to  date 
has  restricted  their  output  power  to  below  IkW  (without  the  use  of  line-broadening 
mechanisms)  due  to  their  limited  mode  size  and  thermal  issues.  Increasing  the  core 
diameter  of  large-mode  area  (LMA)  fibers  beyond  the  commercially  accepted  20-25  pm 
allows  this  limitation  to  be  breached  both  by  using  a  larger  mode  diameter  and  longer 
fibers  for  reduced  thermal  load.  Although  such  very-large-mode-area  (VLMA)  fibers  can 
therefore  allow  the  reduction  of  nonlinear  effects,  their  utility  has  to  date  not  provided 
significant  scaling  due  to  (a)  insufficient  mode  filtering,  and  (b)  unacceptable  mode-area 
reduction  due  to  packaging,  i.e.  bend- induced  mode  deformation  [1,2]. 

Even  more  critical  to  the  degraded  performance  is  the  fiber’s  response  to  thermal 
effects.  First,  the  thermal  gradient  caused  by  pump  absorption  leads  to  a  thermal  lens  that 
inevitably  results  in  spatial  compression  of  the  fundamental  optical  mode  in  the  fiber  [3]. 
Second,  the  dynamics  of  heat  flow  across  the  fiber  cross-section  leads  to  coupling  of 
power  out  of  the  fundamental  mode  into  higher  order  modes  (HOMs),  resulting  in 
significant  dynamic  degradation  of  beam  quality  at  the  output  of  the  fiber  [4]. 

To  date,  comprehensive  modeling  of  large-mode-area  fibers  has  not  yet  been 
addressed  due  to  the  complexity  of  the  problem.  While  it  can  be  envisioned  that  such  a 
comprehensive  model  can  be  computationally  intensive,  the  core  challenges  lie  with  the 
simultaneous  inclusion  of  the  relevant  physical  mechanisms  at  play  in  the  fiber  amplifier: 
multiple  mode  beam  propagation  in  the  fiber,  spatially  resolved  gain  saturation,  a  self- 
consistent  3D  thermal  profile  along  the  cross  section  and  length  of  the  fiber,  and  the 
dynamic  effects  of  mode  instabilities  resulting  from  the  interplay  between  all  of  these 
effects. 

The  goal  of  this  program  is  (a)  to  develop  a  comprehensive  model  for  scaling  the 
mode  area  of  fiber  amplifiers  while  inclusively  addressing  all  of  the  relevant  spatial  and 
dynamic  effects,  and  (b)  to  apply  the  model  to  investigate  VLMA  fiber  amplifiers  scaling 
to  an  effective  SOpm  core  diameter.  The  model  must  produce  results  with  sufficient 
transparency  to  accurately  assess  the  feasibility  of  the  use  of  VLMA  fibers  for  high- 
power  fiber  amplifier  systems. 

2.  Methods,  Assumptions,  and  Procedures 

Two  models  were  specifically  addressed  in  this  program:  a  modally  resolved  model  and  a 
wave  propagation  model.  The  basic  modally  resolved  model  separates  the  power 
evolution  of  a  given  mode  from  its  spatial  field  distribution  [5].  Spatially  dependent  gain 
saturation  is  carried  through  the  population  inversion,  which  is  completely  spatially 
resolved.  Such  a  model  had  already  been  applied  to  assess  beam  quality  in  VLMA 
fibers  [6],  but  required  substantial  development  for  use  in  this  program,  including:  time- 
dependence,  self-consistent  heat  generation,  and  a  (perturbative)  mode-coupling 
mechanism  via  the  thermally  generated  profile.  Although  such  a  model  may  be  lacking 
if  all  of  the  pertinent  physics  are  not  included,  the  advantages  of  this  modally  resolved 
model  are  that  the  physics  is  relatively  transparent  and  the  problem  can  in  principle  be 
analyzed  over  a  very  wide  parameter  space.  More  details  on  this  approach  are  given  in 
Section  4. 
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The  wave  propagation  model  uses  a  three-dimensional  computational  approach. 
This  so-called  “beam  propagation  method”  (BPM)  [7]  does  not  rely  on  modal  resolution 
of  the  optical  beam,  but  rather  calculates  the  evolution  of  the  complete  complex  electric 
field.  The  primary  advantage  of  this  method  is  its  ability  to  model  arbitrary  spatially 
dependent  structures,  such  as  refractive  index  (engineered  or  thermal),  saturable  gain, 
bend  loss,  material  inhomogeneities,  and  a  host  of  other  effects,  while  at  the  same  time 
allowing  for  arbitrary  launch  conditions.  The  output  of  the  model  describes  the  complete 
spatial  profile  of  the  (complex)  optical  field  without  having  a  priori  knowledge  of  the 
modes  of  the  fiber.  The  model  numerically  solves  the  paraxial  wave  equation  with 
transparent  boundary  conditions  to  eliminate  any  real  or  induced  reflections  at  the 
interface  boundaries. 

The  BPM  methodology  had  already  been  successfully  applied  to  VLMA 
fibers  [8],  but  required  two  changes  in  order  to  include  thermal  lensing  and  mode 
instabilities.  The  former  was  accounted  for  by  changing  the  3D  refractive  index  profile 
according  to  the  prescribed  thermal  lens.  The  inclusion  of  mode  instabilities,  however, 
required  that  the  BPM  effectively  become  time  dependent.  Computationally,  this  is 
simply  not  feasible  for  VLMA  fibers.  First,  the  spatial  grid  required  for  computation  is 
already  massive:  200|a,m  x  200)j,m  x  5m.  Second,  the  relative  times  of  interest  are  orders 
of  magnitude  different.  The  optical  time-of-flight  down  a  5m  length  of  fiber  is  about 
25ns.  However,  modal  instabilities  operate  on  a  much  shorter  (~ms)  time  scale.  As  such, 
comprehensive  full-scale  modeling  would  have  required  well  over  10,000  full- fiber 
calculations  before  any  real  changes  become  evident.  Even  using  the  most  sophisticated 
supercomputers,  such  a  calculation  would  take  many  days  to  run,  prohibiting  any 
parametric  analysis  but  the  most  rudimentary  observations.  Instead,  a  different  path  was 
explored  in  which  the  “steady- state”  solutions  were  to  be  applied  iteratively  with  discrete 
jumps  in  time  that  are  significant  compared  to  the  time-of-flight  down  the  fiber,  but 
minute  compared  to  the  time  scale  of  the  dynamical  instabilities  under  study.  In  this  way, 
the  full-scale  BPM  strategy  could  be  applied  to  the  dynamic  problem  while 
simultaneously  minimizing  computation  time  to  a  reasonable  level  suitable  for  studying 
the  problem  at  hand.  More  details  on  this  approach  are  given  in  Section  3. 

Finally,  it  was  noted  throughout  the  program  that  even  using  the  Blue  Hive 
supercomputer  at  the  University  of  Rochester,  the  simulations  took  exceedingly  long 
periods  of  time  to  execute.  For  example,  using  our  modified  modal  model  to  find  the 
instability  threshold  for  the  simplest  case  studied  (two  modes  in  a  20pm  LMA  fiber) 
would  take  24  hours  to  complete  using  our  fiilly  allotted  24  nodes.  Therefore,  a 
methodology  was  developed  to  assess  the  viability  of  specific  fiber  designs  before 
running  the  detailed  simulations.  Specifically,  the  data  resulting  from  the  bulk  of  our 
initial  runs  was  used  to  generate  a  set  of  engineering  rules  for  designing  and  down- 
selecting  specific  fiber  types  without  using  any  propagation  or  thermal  modeling.  These 
engineering  rules  allowed  us  to  eliminate  non-viable  solutions  without  wasting  valuable 
computational  resources,  and  focus  on  the  development  of  promising  new  VLMA  fiber 
types.  More  details  will  be  presented  in  Section  6. 

It  must  be  noted  that  throughout  our  analysis,  the  fiber  types  analyzed  were 
focused  on  conventional-type  fibers  (including  polarization  maintaining  fibers)  using 
waveguide-based  physics,  i.e.  no  photonic  crystal  or  photonic  band-gap  fiber  types  were 
investigated.  Ultimately,  such  alternative  fiber  types  become  challenging  to  fabricate  to 
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spec  and  are  highly  susceptible  to  packaging  (i.e.,  bending)  and  external  thermal 
influences. 


3.  Beam  Propagation  Model  Development 

One  of  the  key  aspects  of  the  BPM  model  is  its  ability  to  indiscriminately  propagate  the 
entire  complex  optical  field  rather  than  relying  on  modal  decompositions.  We  define  this 
aspect  as  critical  since  the  large-scale  temperature  changes  inside  the  pumped  fiber 
dictate  that  the  cold-fiber  modes  are  insufficient  to  describe  the  propagation  of  light  down 
the  fiber  from  a  strictly  perturbative  perspective,  as  is  commonly  done  in  modally 
resolved  models.  This  is  particularly  relevant  since  the  temperature  profile  is  not  even 
constant  along  the  length  of  the  fiber,  especially  in  the  case  of  single-end  pumping. 

3.1  Full  BPM  Model 

Our  BPM  model  numerically  solves  the  paraxial  wave  equation  for  the  optical  (electric) 
field  E  along  the  propagation  (z)  axis  as 
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where  E  is  the  optical  field,  where  ko  =  1%I\  and  nR  is  a  reference  refractive  index.  On 
the  right  hand  side  of  Eq.  (1),  the  first  term  represents  diffraction,  and  the  second  term 
accounts  for  the  refractive  index  profile,  n(x,y,z),  in  three  dimensions,  including  the 
engineered  cross-sectional  fiber  profile  (constant  in  z)  and  the  thermal  contribution.  The 
last  term  incorporates  the  spatial  dependence  of  the  gain  doping  profile  through  go(x,y) 
and  accounts  for  spatially  localized  saturable  gain  in  three  dimensions,  fat  is  the 
saturation  intensity  obtained  including  the  influence  of  pump  saturation  [9]  through  a 
standard  ytterbium  (Yb)  kinetics  model  [10]. 

Eq.  (1)  is  solved  with  a  finite-difference  scheme  (FD-BPM)  [11]  using  the 
Alternating  Direction  Implicit  (ADI)  method  [12].  The  need  to  model  long  propagation 
lengths  required  that  transparent  boundary  conditions  be  implemented  to  eliminate  any 
real  or  induced  reflections  at  the  interface  boundaries  [13]. 

As  stated  earlier,  the  methodology  for  applying  the  BPM  to  the  dynamic  thermal 
model  is  to  iteratively  calculate  the  optical  and  thermal  properties  of  the  amplifier.  First, 
the  speed-of-light  (ns-scale)  propagation  of  the  optical  beam  is  calculated  through  the 
fiber  using  Eq.  1 .  Then,  the  full  3D  thermal  profile  in  the  fiber  is  calculated  using  the 
standard  2D  time-dependent  heat  equation  at  each  cross-sectional  step  along  the  fiber 
using  heating  that  is  spatially  defined  via  the  pump/signal  energy  deposition/extraction. 
The  dimensional  reduction  (3D  to  2D)  is  reasonably  valid  due  to  the  adiabatic  nature  of 
the  average  thermal  evolution  along  the  fiber,  and  reduces  computation  time 
significantly.  This  two-step  process  (optical  /  thermal)  is  repeated  on  thermal  (ps) 
timescales  in  order  to  simulate  the  dynamics  represented  by  the  thermal  mode  instability. 

Reasonably  early  in  the  program,  the  results  indicated  that  the  computation  is  too 
memory  intensive  and  insufficiently  fast  for  anything  but  a  few  singular  demonstrative 
simulations.  The  computation  grid  was  approximately  200  x  200  x  250,000  =  10  billion 
data  points  each  in  optical  amplitude,  optical  phase,  and  temperature  (which  must  be 
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retained  from  one  iteration  to  the  next).  Although  it  was  initially  suspeeted  that  the 
temperature  grid  could  be  reduced  by  a  significant  fraction  (~100x),  it  was  recognized 
that  this  would  only  be  true  if  one  considers  only  the  two  lowest-order  modes  for  a  low 
numerical  aperture  (NA  =  0.06)  geometry.  Since  some  of  our  fiber  concepts  indicated 
that  the  second  (LPn)  mode  is  effectively  nonexistent  and  others  require  the  use  of  higher 
NA  cores,  neither  of  the  requirements  above  are  met  universally,  and  the  entire  grid  must 
therefore  be  calculated  (and  retained)  in  general.  As  such,  all  data  must  be  kept  at  the  3  x 
lOB  data  size,  which  is  simply  prohibitive. 

3.2  Reduced  BPM  Model 

Two  methods  for  minimizing  the  computation  time  and  data  volume  were  investigated. 
The  first  method  was  reducing  the  spatial  grid  (both  transverse  and  axial  since  they  are 
effectively  coupled  in  order  to  produce  accurate  BPM  results).  If,  for  example,  a  factor 
of  5  can  be  gained  in  each  of  the  three  dimensions,  then  the  three  arrays  (optical 
amplitude,  optical  phase,  and  temperature)  can  be  reduced  to  SOM  data  points  each, 
which  would  represent  a  viable  data  volume  for  the  required  exchange.  Using  this 
approach,  however,  appropriate  care  must  be  taken  to  ensure  preservation  of  all  beat 
lengths  for  proper  validity  of  the  model. 

In  the  cases  that  were  modeled  using  reduced-grid  conditions,  it  was  discovered 
that  although  reducing  the  resolution  sufficiently  to  allow  reasonable  data  sizes  and 
computation  time  resulted  in  data  that  reproduced  similar  physical  behavior  under  certain 
conditions,  it  was  sufficiently  inaccurate  for  use  as  a  predictive  tool  for  designing  and 
engineering  new  fiber  types,  even  using  focusing  on  the  LPoi  and  LPn  modes.  As  a 
simple  example,  the  amplifier  efficiency  should  not  change  considerably  under 
dynamically  instable  conditions  compared  to  stable,  continuous-wave  (CW)  operation. 
However,  in  some  cases  that  were  modeled  with  a  meaningfully  reduced  spatial  grid,  the 
amplifier  efficiency  was  too  high  to  be  physically  realistic.  Due  to  this  and  other 
observed  inconsistencies  and  inaccuracies,  reducing  the  spatial  grid  resolution  was 
deemed  as  a  non- viable  option  for  minimizing  data  volume  and  computation  time. 

The  next  path  attempted  was  to  implement  a  simplified  thermal  predictor  model 
that  would  not  require  the  storage  and  transfer  of  huge  data  volumes.  In  general,  thermal 
effects  are  included  in  the  model  via  induced  change  to  the  3D  refractive  index  profile 
according  to  the  differential  form  An(x,y,z)  =  (dn/dT)AT(x,y,z),  where  dn/dT  is  the 
change  of  the  refractive  index  of  glass  with  respect  to  temperature  (~10“^/°C).  The 
localized  temperature  change  AT(x,y,z)  is  essentially  driven  by  energy  exchange  between 
the  pump  and  signal  photons  within  the  medium,  and  is  conventionally  calculated  via  the 
heat  equation.  The  heat  load  that  drives  these  thermal  changes  is  given  by 
(l-'nQD)Ps(OeN2-aaNi),  whcrc  (l-riQo)  is  the  quantum  defect,  Ps(x,y,z)  is  the  spatially 
resolved  signal  power,  Oe/a  is  the  emission/absorption  cross  section  at  the  signal 
wavelength,  and  N2/i(x,y,z)  is  the  spatially  resolved  upper/lower  state  electron  density. 
Note  that  the  non-radiative  decay  (N2/T2)  has  been  neglected  for  simplicity  since  its 
impact  should  be  negligible  in  a  high-efficiency  amplifier. 

This  formalism  can  be  further  simplified  by  recognizing  that  (aeN2-aaNi) 
represents  the  localized  gain  profile.  In  this  case,  the  heat  load  can  then  be  simplified  as 
(l-'nQD)Ps(x,y,z)g(x,y,z),  where  the  spatial  dependence  has  been  explicitly  shown.  Since 
the  optical  profile  is  directly  linked  to  the  heat  load  through  the  equation  described  above, 
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this  means  that  the  heat  load  will  change  before  the  heat  can  be  redistributed.  As  a  coarse 
approximation,  this  means  the  temperature  profile  may  be  considered  to  be  nearly 
equivalent  to  the  heat  load  for  the  purpose  of  implementing  our  serial  BPM  scheme. 

Eliminating  the  calculation  of  the  3D  temperature  profile  along  the  entire  can 
result  in  significant  savings  in  terms  of  computation  time.  However,  the  translation  from 
heat  load  to  temperature  is  a  critical  step  that  is  omitted  in  the  simplified  formalism 
outlined  above.  Since  conductive  heat  transfer  is  being  ignored,  straightforward 
conversion  from  heat  to  temperature  will  result  in  higher  temperatures  than  expected,  and 
the  ps-scale  diffusion  implicit  in  the  heat  equation  solutions  will  be  absent. 
Implementing  this  approximation  will  therefore  have  an  exaggerated  impact  on  the 
optical  field,  which  will  then  generate  a  substantially  different  heat  load,  resulting  in  a 
rapidly  (and  non-physically)  varying  thermal  profile. 

To  mitigate  this  problem  and  emulate  thermal  diffusion,  the  process  that  was 
implemented  consisted  of  mixing  the  newly  generated  thermal  profile  with  the  prior 
iteration  thermal  profile  using  a  relative  strength  of  a  <  1.  The  original  idea  was  to 
determine  the  value  of  a  via  merit  function  using  an  adiabatic  requirement  on  the 
evolution  of  the  thermal  profile.  However,  we  discovered  through  our  simulations  that 
adiabatic  evolution  of  the  temperature  profile  required  a  value  of  a  that  was  simply  too 
low  to  generate  mode  instabilities.  Instead,  omitting  the  adiabatic  requirement  and  using 
larger  values  of  a  resulted  in  dynamic  spatial  instabilities.  However,  the  thermal  profiles 
in  this  case  changed  in  an  unphysical  manner  and  were  thus  deemed  unreliable  from  a 
physics  perspective. 

3.3  Final  BPM  Model  Assessment 

The  result  of  our  efforts  with  the  BPM  modeling  is  that  the  method  is  either  (a)  too  large 
and  computationally  intensive  to  represent  anything  other  than  a  singular  demonstration 
(which  has  already  been  published  by  other  groups  [14]),  or  (b)  too  inaccurate  to  produce 
meaningful  trends  or  design  validation.  These  findings  therefore  eliminate  iterative  BPM 
as  a  viable  design  tool  for  including  the  prediction  of  dynamic  thermal  mode  instabilities. 

4.  Rate-Equation  Model  Development 

The  basic  modally  resolved  model  separates  the  power  and  phase  evolution  of  a  given 
mode  from  its  spatial  field  distribution.  One  of  the  fundamental  drawbacks  of  this 
method  is  that  the  modes  are  fixed  through  the  modeling  by  de  facto,  with  only 
perturbative  coupling  between  the  modes  to  modify  their  behavior  (power/phase),  in 
order  to  facilitate  rapid  simulation.  For  the  case  of  asymmetrically  pumped  high-power 
fiber  amplifiers,  this  represents  a  flaw  in  using  a  modally  decomposed  model,  and  was 
our  original  motivation  for  attempting  to  apply  the  BPM  methodology  to  this  problem. 
However,  we  have  developed  a  new  method  that  overcomes  this  fundamental  limitation 
while  still  retaining  all  of  the  benefits  of  modal  decomposition  modeling,  as  will  be 
described  in  Section  4.2. 

The  rate  equation  modeling  (REM)  work  started  with  an  analysis  of  all  relevant 
existing  models.  A  dynamic  model  in  particular  was  selected  since  non-dynamic  models 
are  incomplete  by  their  very  nature.  While  steady-state  and  stability-analysis  type 
models  are  typically  excellent  at  predicting  trends,  they  are  also  notorious  for  requiring 
arbitrary  fit  parameters  to  match  experimental  data.  Therefore,  only  dynamic  models 
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were  considered.  Although  many  REM  models  exist  (e.g.,  DTU,  Jena,  Clemson,  AS 
Photonics,  U.  Mich.,  etc.),  the  model  from  the  Air  Force  Research  Laboratory 
(AFRL)  [15]  was  down-selected  as  our  initial  starting  point.  AFRL’s  modally 
decomposed  REM  not  only  includes  all  of  the  relevant  physics  in  a  true  dynamic  model, 
but  it  does  so  in  a  transparent  (physics-ready)  methodology  derived  from  fundamental 
(Maxwell’s)  equations.  Perhaps  most  importantly,  it  allows  modification  for  inclusion  of 
higher-order  physics.  It  is  also  very  similar  to  our  existing  steady-state  modally 
decomposed  REM  model  and  allows  for  simpler  coding  transition. 

4.1  Basic  REM  Model  Definition  and  Validation 

Modal  propagation  down  the  fiber  amplifier,  incorporating  coupling  between  modes  due 
to  refractive  index  perturbations  and  amplifier  gain,  is  described  by  the  singular  and 
deceptively  simple  equation 


dz 


(2) 


where  Ei  is  the  complex  field  of  the  i*  mode,  and  (3i  is  its  propagation  coefficient.  Apart 
from  the  summation  over  all  modes,  Eq.  2  looks  similar  to  a  single-mode  propagation 
equation  where  the  kappas  take  the  place  of  index  and  gain  respectively.  Their 
mathematical  forms  are  given  by 

(^)  =  (X’y)Sn(x,y,z)a)^  (x,y)dxdy  (3) 

’^j,k(z)  =  JJ  ^i(x,y)[c^sN2(x,y.z)-c?sNi(x,y,z)]Oi^(x,y)dxdy  (4) 


where  <l>k(x,y)  is  the  power-normalized  spatial  mode  profile,  6n(x,y,z)  is  the  modification 
to  the  base  refractive  index  (embedded  in  p),  are  the  emission/absorption  cross- 
sections,  N2/i(x,y,z)  are  the  excited/ground  state  population  inversion  densities,  related  to 
the  total  doping  density  via  NT(x,y)  =  N2(x,y,z)  +  Ni(x,y,z). 

For  the  case  of  only  a  single  mode,  Eq.  3  represents  a  modification  to  the 
propagation  coefficient  Pj,  Eq.  4  represents  gain,  and  Eq.  2  reduces  to  the  familiar 
propagation  equation.  In  the  multimode  case,  this  still  holds  true  for  the  j  =  k  terms, 
however  there  exist  cross  terms  that  determine  the  coupling  between  modes.  As  in 
conventional  waveguide  theory,  this  coupling  is  determined  by  the  spatial  overlap  of  the 
two  modes  over  the  region  of  the  perturbation,  in  this  case  the  modification  to  the 
refractive  index  and  the  material  gain.  Note  that  the  gain  and  index  change  over  the 
length  of  the  fiber,  resulting  in  coupling  coefficients  that  are  z-dependent. 

The  time-dependent  heat  equation  is  solved  between  iterations  of  the  optical 
solutions  to  determine  the  refractive  index  perturbation  Sn(x,y,z)  at  each  point  along  the 
fiber.  The  populations  N2/i(x,y,z)  are  determined  by  spatially  dependent  ytterbium 
kinetics  rate  equations  using  the  spatial  distributions  of  the  signal  and  pump  power  at 
each  point  along  the  fiber.  The  pump  power  is  assumed  to  be  spatially  uniform  and  is 
calculated  via  the  spatially  integrated  populations. 

Although  the  model  does  not  contain  noise  to  initiate  the  instability,  it  is  in  fact 
not  required.  No  practical  seed  launch  into  the  amplifier  is  100%  efficient  in  launching 
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into  the  fundamental  mode  of  the  fiber,  and  power  will  be  injeeted  (unintentionally)  into 
other  modes  of  the  fiber.  Even  if  small,  the  seed  power  injeeted  into  other  modes  will 
dominate  quantum  noise  (spontaneous  emission),  thus  obviating  the  need  for  inclusion  of 
noise  in  the  model. 

The  model  was  validated  by  simulating  fiber  amplifiers  two  meters  in  length  with 
a  20pm  core  diameter  and  a  230pm  pump  cladding.  The  nominal  (room  temperature) 
cladding  refractive  index  was  taken  to  be  1.45  with  a  core  NA  of  0.06.  The  doping 
density  of  the  core  was  taken  as  6.2  xlO  m  ,  with  pump  and  signal  wavelengths  of 
975  nm  and  1080  nm,  respectively.  Other  parameters  were  nominally  selected  for  a 
convectively  cooled  Yb-doped  silica  fiber.  The  axial/temporal  step  size  was  determined 
to  sufficiently  resolve  the  mode  beating  while  simultaneously  ensuring  conservation  of 
energy  in  the  system.  Although  smaller  step  sizes  were  also  investigated,  the  results  were 
similar  but  at  the  expense  of  additional  computation  time. 

Though  the  model  describes  the  coupling  between  all  supported  fiber  modes, 
coupling  between  the  fundamental  and  the  next  highest  order  mode  is  the  strongest  and 
thus  will  determine  the  thermal  most  instability  (TMI)  threshold  in  standard  (step-index) 
fiber  types.  For  this  kW-class  amplifier,  35  W  of  signal  power  was  seeded  into  the 
fundamental  (LPoi)  mode,  with  1.75  W  (5%  of  the  fundamental-mode  power)  launched 
into  the  higher-order  (LPn)  mode.  The  pump  was  co-propagated  with  the  signal  launch. 

Figure  1  shows  the  power  evolution  along  the  length  of  the  FMA  fiber  amplifier 
of  the  FPoi  and  FPn  modes  at  sequential  steps  in  time  for  several  pump  powers.  At 
750  W  pump  power,  the  amplifier  shows  stable  behavior  at  all  times  (noted  by  the  fact 
that  the  FPoi  and  FPn  curves  are  all  identical  for  each  time  step,  as  indicated  in  the  plot), 
with  the  FPoi  mode  containing  most  of  the  optical  power.  At  850  W  pump  power,  modal 
power  interchange  along  the  fiber  occurs  and  is  dynamic,  noted  by  the  fact  that  the  modal 
curves  representing  the  various  time  steps  are  not  identical.  At  950  W  pump  power,  the 
modal  power  transfer  initiates  earlier  along  the  fiber,  indicating  a  stronger  instability. 
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Figure  1:  LPgi  and  LPn  power  evolution  along  the  fiber  amplifier  length  at  various  steps  in  time. 
The  seed  power  is  35  W  and  the  pump  power  is  show  in  each  figure  inset.  Various  time  steps  are 
shown  as  additional  curves  in  the  plots,  with  the  t  =  0  case  shown  as  green  and  yellow  for  the 
LPoi  and  LPn  modes  respectively. 


Figure  2  shows  the  temporal  evolution  of  the  power  in  each  mode  (LPoi  and  LPn) 
at  various  pump  powers.  Each  simulation  starts  with  a  dynamic  power  transfer  between 
the  modes  that,  when  under  the  instability  threshold,  stabilizes  to  CW-type  conditions 
within  30  ms;  this  represents  the  settling  time  of  the  simulation  and  is  dependent  on  the 
initial  thermal  profile  selected.  As  the  pump  power  is  increased,  the  results  indicate  the 
tendency  towards  instability  as  the  settling  time  (to  CW-type  conditions)  becomes  larger 
with  increasing  pump  (and  thus  output)  power.  Dynamic  instability  is  finally  reached  in 
the  range  750-800W  pump  power,  beyond  which  the  power  exchange  does  not  stabilize. 

The  minute  fluctuations  near  threshold  (see  the  blue  curve  at  750W  pump)  are  key 
for  validating  the  (temporal)  length  of  runs  required,  as  the  rapid  power  transfer  between 
modes  ceases  below  the  instability  threshold,  but  can  be  observed  as  threshold  is 
approached. 
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Figure  2:  LPqi  (blue)  and  LPn  (orange)  power  at  the  amplifier  output  as  a  function  of  time  over  a 
35-ms  time  span.  The  seed  power  is  35  W  and  the  pump  power  is  inset  in  each  plot. 


Similar  runs  were  executed  for  various  seed  powers,  resulting  in  the  (coarse) 
instability  thresholds  shown  in  Fig.  3.  The  data  shows  that  the  instability  threshold 
increases  with  seed  power,  which  is  not  only  physically  justifiable  (due  to  the  increased 
saturation  by  the  fundamental  mode)  but  was  also  previously  predicted  [15]. 
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Figure  3:  Threshold  pump  power  for  the  onset  of  thermally  induced  mode  instability  as  a  function 
of  seed  power. 

The  results  of  Figures  1-3  demonstrate  that  the  base  model  as  coded  yields 
predictions  that  reproduce  the  requisite  physics  at  a  level  that  is  sufficiently  quantitative 
for  fiber  design  and  engineering  of  realistic  high-power  fiber  amplifiers. 

4.2  Modification  to  the  REM  Model 

As  stated  previously,  the  fundamental  flaw  of  using  modally  resolved  models  for  our 
application  is  that  the  mode  and  its  eigenvalue  ((3)  are  fixed  throughout  the  entire 
simulation  space  (both  in  time  and  along  the  length  of  the  fiber).  However,  the  varying 
thermal  characteristic  along  the  length  of  the  fiber,  given  by  Eqs.  (2)-(4),  means  that  the 
natural  mode  of  the  fiber  (i.e.,  the  ideal  waveguide  mode  solution)  can  also  change  as  a 
function  of  the  length  along  the  fiber  and  in  time  as  the  thermal  characteristic 
dynamically  changes.  Since  the  pump  power  varies  along  the  fiber  length,  drastically 
different  thermal  profiles  from  one  end  of  the  fiber  to  the  other  will  lead  to  different 
transverse  index  profiles  along  the  fiber.  This  directly  translates  modes  having  different 
transverse  profiles  on  opposite  ends  of  the  fiber. 

Provided  the  thermal  changes  are  sufficiently  slow  compared  to  the  mode 
formation  length,  any  given  mode  of  the  fiber  will  adiabatically  change  along  the  length 
of  the  fiber  with  minimal  coupling  between  modes  (solely  due  to  the  axial  temperature 
changes,  and  not  including  transverse  effects).  It  is  precisely  this  effect  that  we  add  into 
our  REM  model.  Specifically,  Oj(x,y)  becomes  Oj(x,y,z),  and  (3j  becomes  pj(z).  Note 
that  this  not  only  impacts  the  coupling  coefficients  in  Eqs.  (3)  and  (4),  but  also  the  phase 
matching  condition  in  Eq.  2. 

In  this  new  formalism,  we  calculate  the  modes  of  a  given  fiber  as  a  function  of 
temperature,  which  manifests  in  the  well-known  graded-index  profile  for  cylindrical 
symmetries  [16].  The  sets  of  modes  (one  set  of  unique  transverse  modes  at  each 
temperature)  are  solved  using  an  existing  finite-element-analysis  (FEA)  mode-solver.  In 
this  way,  any  given  mode  (e.g.,  LPoi,  LPn,  etc.)  has  both  a  spatial  profile  and  a 
propagation  coefficient  that  vary  (discretely)  as  a  function  of  core  temperature. 
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These  results  are  applied  to  the  REM  model  as  follows:  onee  the  spatial 
temperature  profile  is  ealculated  at  a  given  point  in  the  fiber,  the  resulting  average  eore 
temperature  is  used  to  interpolate  (a)  the  shape  of  eaeh  mode  at  that  point  in  the  fiber 
(using  the  previously  solved  mode  set)  and  (b)  their  propagation  eoefficients  pi.  Akin  to 
the  adiabatic  physics  in  the  actual  (physical)  fiber  amplifier,  adiabatic  conversion  with 
propagation  is  assumed  such  that  the  power  does  not  transfer  between  modes;  rather  the 
mode’s  spatial  profile  and  effective  index  change  as  the  mode  propagates  down  the  fiber 
without  transfer  of  power  to  any  other  modes  (neglecting  transverse  effects).  In  this  way, 
the  evolution  of  the  mode  shapes  and  the  beating  between  the  modes  are  taken  into 
account  in  a  natural  and  physical  way,  and  the  remaining  (transverse)  thermal 
perturbations  can  be  truly  perturbative. 

Figure  4  shows  lineouts  of  the  LPoi  and  LPn  modes  calculated  for  both  a  standard 
25|am,  0.06-NA  LMA  fiber  and  a  similar  fiber  with  SOpm  core  diameter.  Each  plot 
shows  the  overlap  of  many  modes  at  different  (peak)  core  temperatures  ranging  from  0  to 
900°C  above  ambient.  For  this  relatively  small  step-index  fiber,  the  mode  shape  is  only 
very  weakly  impacted  (compressed)  by  the  approximately  parabolic  temperature  profile 
rising  in  the  core  at  high  powers. 


(a) 


Transverse  Coordinate  (|jni) 
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Figure  4:  LPqi  and  LPn  mode  lineouts  for  0.06  NA  fiber  cores  with  diameter  of  (a)  25  jum  and 
(b)  50  jim.  Multiple  curves  on  each  plot  range  from  0  (blue)  to  900°C  (orange)  peak  temperature 
above  ambient. 


Figure  5  shows  the  difference  in  the  effective  index  between  the  LPoi  and  LPn 
propagation  modes  as  a  function  of  temperature  for  the  same  25 pm  and  50pm  fibers.  It 
is  important  to  note  that  for  both  fibers,  the  increased  thermal  load  actually  causes  an 
increase  in  the  separation  of  the  mode  eigenvalues,  which  would  nominally  lead  to  a 
decrease  in  the  phase  matching  [see  Eq.  (2)].  Although  the  relative  increases  are  small 
(2.2%  for  the  25pm  core,  and  5.4%  for  the  50pm  core),  the  additional  complications  in 
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the  simulation  make  the  results  more  reliable  for  use  in  the  development  of  new  VLMA 
fibers. 


Temperature  rise  above  ambient  (®C) 


Figure  5:  Effective  index  difference  between  LPqi  and  LPn  modes  for  0.06  NA  fiber  cores  with 
diameter  of  (a)  25  jum  and  (b)  50  jum. 

The  inelusion  of  the  thermally  varying  mode  profiles  and  propagation  eoefficients 
were  subsequently  used  in  all  cases  moving  forward  in  order  to  provide  a  more  accurate 
quantitative  description  of  the  physics  resulting  from  the  various  fiber  designs  under 
consideration  in  the  VLMA  analysis. 
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5.  Core  Scaling  with  Conventional  Step-Index  Fibers 

The  REM  model  described  in  Sections  4  were  applied  to  various  fiber  amplifiers  to 
obtain  a  deeper  understanding  of  the  physics  and  understand  the  impact  of  thermal  mode 
instability  on  core-area  scaling.  Unless  otherwise  stated,  all  parameters  are  the  same  as 
described  in  Section  4. 

One  interesting  aspect  allowed  by  the  REM  simulations  is  the  ability  to  observe 
the  transition  between  stable  and  unstable  operation.  Figure  6  shows  a  comparison  of 
longitudinal  temperature  profiles  taken  7  pm  off  center  of  the  fiber  propagation  axis  over 
the  last  8  mm  of  the  fiber.  The  below-threshold  curve  shows  a  periodic-type  temperature 
variation,  indicating  interferometric  beating  (without  power  exchange)  between  the  high- 
power  LPoi  and  low-power  LPn  modes.  In  contrast,  the  above-threshold  curve  shows  a 
much  more  complicated  temperature  profile,  indicating  significant  modal  coupling. 


Distance  from  end  of  liber  (mm) 

Figure  6:  Temperature  variation  above  the  local  temperature  7  pm  off  the  center  of  the  fiber  axis. 
Blue  and  red  curves  are  below  and  above  threshold,  respectively. 

5.1  Core-Diameter  Scaling 

The  nature  of  TMI  is  insidious  in  that  even  though  beam  quality  can  be  retained  at  low 
power,  LMA  fibers  driven  to  higher  powers  exhibit  thermally  induced  mode  instability 
[17]  that  is  characterized  by  a  significant  fraction  of  optical  power  being  exchanged 
between  modes  of  the  fiber  at  kHz  rates.  To  date,  some  methods  have  been  demonstrated 
to  reduce  the  impact  of  TMI  and  increase  its  threshold  [18],  which  again  leads  to  the 
exploration  of  fiber  amplifiers  with  larger  core  areas.  It  is  therefore  prudent  to 
understand  what  happens  to  the  TMI  threshold  when  the  core  area  is  increased  for 
conventional  step-index  fiber.  Although  some  preliminary  modeling  on  step  index  fibers 
has  been  carried  out  [19],  the  work  used  a  significantly  simplified  model  that  while 
saving  computation  time,  neglected  critical  physics  such  as  transverse  spatial  hole 
burning  which  is  at  the  core  of  TMI  physics.  Moreover,  the  results  are  dependent  on  an 
arbitrary  frequency  difference  between  the  modes  that  not  only  obscures  the  analysis,  but 
is  also  unphysical  from  the  perspective  of  realistic  launch  conditions. 
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The  REM  model  was  applied  to  fibers  with  core  diameters  varying  from  20  to 
50  pm,  each  with  230pm  cladding  diameter.  The  doping  density  of  the  core  was  adjusted 
between  2.5  and  6.2  x  10  m  to  maintain  the  same  fiber  length  and  output  power  for  a 
given  pump  power  as  the  core  diameter  is  adjusted.  Simulations  were  run  with  25  W 
increments  in  pump  power,  and  the  TMI  threshold  was  determined  when  dynamic  power 
exchange  between  modes  was  observed  (as  shown  in  Fig.  2). 

Figure  7  shows  the  resulting  TMI  threshold  for  core  diameters  ranging  from  20  pm  to 
50  pm  for  forward  a  forward-pumped  amplifier.  The  error  bars  indicate  the  25  W  pump 
power  increment.  Figure  8  shows  a  monotonically  decreasing  TMI  threshold  with 
increasing  core  diameter.  At  first  glance,  this  may  seem  counter  intuitive  since  the  heat  is 
generated  over  a  larger  area.  However,  thermal  effects  (specifically  thermal  lensing) 
have  been  observed  to  have  larger  impact  on  larger  core  diameters  [20],  a  fact  that  has 
also  been  confirmed  via  analytic  and  numerical  modeling  [21]. 
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Figure  1:  TMI  threshold  for  forward  pumped  LMA  fiber  as  a  function  of  core  diameter.  The  line 
is  a  guide  for  the  eye,  and  the  error  bars  represent  the  25  W  pump  increment. 


The  trend  observed  in  Fig.  7  means  that  even  though  TMI  currently  represents  a 
challenging  in  power  scaling  LMA  fiber  amplifiers,  it  also  represents  a  significant  barrier 
to  increasing  the  core  diameter  for  mitigating  nonlinear  effects  unless  other 
methodologies  are  brought  in  to  combat  TMI.  This  will  be  the  focus  of  the  remainder  of 
this  work.  As  such,  Fig.  7  is  not  only  valuable  in  physical  understanding  of  the 
phenomenon,  but  also  in  comparing  new  fiber  designs  to  conventional  step-index  fiber. 

5.2  Pumping  Considerations 

Since  the  thermal  issues  resulting  in  TMI  are  due  to  heat  distributed  by  pump  absorption, 
it  is  natural  to  consider  the  impact  of  pump  geometry  on  the  TMI  threshold.  Figure  8 
represents  the  typical  case  of  forward  pumping  (i.e.,  co-propagation  the  pump  and  signal 
light).  Notionally,  this  places  a  large  heat  burden  towards  the  input  end  of  the  fiber.  The 
opposite  case,  backward  pumping  (i.e.,  counter-propagating  the  pump  and  signal  light) 
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places  a  large  heat  burden  towards  the  output  end  of  the  fiber.  Therefore,  one  may 
consider  that  backward  pumping  will  not  result  in  any  improvement  (i.e.,  increase  in  TMI 
threshold).  However,  if  we  instead  consider  distributing  the  heat  burden  throughout  the 
fiber,  for  example  by  including  both  forward  and  backward  pumping,  then  it  may  be 
reasonable  to  expect  an  improvement  to  manifest  as  increase  in  the  TMI  threshold. 

Figure  8  shows  the  temperature  distribution  as  a  function  of  length  along  the  fiber 
for  three  cases:  forward  pumped,  backward  pumped,  and  bi-directionally  pumped.  The 
bi-directional  pumping  scheme  was  specifically  designed  to  flatten  out  the  temperature 
curve  as  much  as  possible  with  the  goal  of  increasing  the  TMI  threshold. 


Figure  8:  Thermal  distribution  along  25-jum  LMA  fiber  for  forward  (blue),  backward  (orange), 
and  bi-directional  (purple)  pumping. 

Figure  9  shows  the  LPoi  power  along  the  length  of  the  fiber  at  the  TMI  threshold 
for  each  case.  As  might  be  expected  at  TMI  threshold,  the  mode  instability  occurs  within 
the  fiber  where  the  temperature  reaches  its  highest  point.  For  the  bi-directionally  pumped 
case,  although  the  temperature  profile  is  relatively  flat,  its  peak  temperature  occurs  only  a 
little  further  along  the  fiber  than  for  the  forward-pumped  case  and  thus  the  modal  power 
exchange  occurs  only  a  little  farther  down  the  fiber  (albeit  at  higher  pump  power). 

The  flatness  of  the  bi-directionally  pumped  case,  however,  does  have  its  benefits. 
One  can  immediately  see  from  Fig.  9  that  the  temperature  peak  in  the  bi-directional  case 
is  lower  than  that  of  the  forward  pumped  case.  It  stands  to  reason  that  the  TMI  threshold 
would  be  correspondingly  higher,  which  was  borne  out  numerically.  Figure  10  shows  the 
results  for  the  three  pumping  cases.  In  the  bi-directionally  pumped  case,  the  fraction  of 
forward/backward  pumping  was  optimized  at  each  core  diameter  to  provide  the  flattest 
thermal  profile  along  the  fiber. 
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Figure  9:  LPqi  mode  power  distribution  along  25pm  LMA  fiber  for  forward  (blue),  backward 
(orange),  and  bi-directional  (purple)  pumping,  at  their  respective  TMI  thresholds. 


Fiber  Core  Diameter  (fim) 

Figure  10:  TMI  threshold  for  forward  (blue),  backward  (red),  and  bi-directional  (green)  pumping 
as  a  function  of  core  diameter.  The  lines  are  guides  for  the  eye,  and  the  error  bars  represent  the 
pump  increment. 

Although  the  backward-pumped  case  does  not  provide  any  benefit  due  to  the  very 
high  temperatures  at  the  end  of  the  fiber,  the  bi-directionally  pumped  case  offers  a  ~10% 
increase  in  the  TMI  threshold.  Even  though  in  itself  this  benefit  is  not  substantial,  it  may 
be  more  beneficial  when  combined  with  other  TMI-mitigation  methods.  Note  that  this 
improvement  is  less  than  that  predicted  using  the  BPM  method  [22],  although  the  initial 
conditions  for  those  simulations  were  not  stated.  Other  researchers  have  investigated 
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double  passing  the  signal,  in  both  amplifier  and  laser  eonfigurations,  and  found  a 
reduetion  in  the  TMl  threshold  in  both  cases  compared  to  that  of  the  single-pass 
amplifier  [23]. 

The  results  of  this  section  make  it  abundantly  clear  that  new  methods  for  scaling 
the  core  area  of  LMA  fibers  need  to  be  developed.  This  will  be  the  focus  of  the 
remainder  of  the  document. 

6.  Novel  VLMA  Fiber  Designs 

The  results  of  the  prior  section  demonstrate  that  is  it  critical  to  seek  new  fiber  designs  for 
core-area  scaling.  Nominally  speaking,  the  analysis  of  VLMA  fibers  can  proceed  straight 
forwardly  using  the  REM  model.  However,  computational  resources  set  practical  limits 
on  what  can  be  achieved  in  a  reasonable  amount  of  time.  Specifically,  a  single  set  of 
two-mode-only  runs  with  ramped  pump  power,  each  with  sufficient  time  to  observe 
modal  instability  and  its  potential  stabilization,  (i.e.,  one  set  of  runs  to  determine  the 
instability  threshold)  requires  ~24  hours  of  run  time  exploiting  120  processors  on  24 
nodes  using  the  University’s  Blue  Hive  supercomputer.  In  addition,  certain  designs 
require  multiple  mode  pairs  (e.g.  LPoi-LPn,  LP01-LP02,  etc.)  to  be  investigated  to 
determine  which  mode  has  the  lowest  instability  threshold.  Therefore,  preliminary 
assessment  methods  have  been  developed  to  minimize  computation  time. 

Although  the  program  was  originally  stated  to  explore  only  conventional  step- 
index  fibers  using  gain  filtering  (confined  gain),  inclusion  of  further  index  profile 
engineering  was  viewed  as  a  significant  opportunity  not  to  be  missed.  Only  conventional 
index-profiling  designs  were  considered;  photonic  crystal  fiber  and  photonic  band-gap 
fiber  were  deemed  outside  the  range  of  interest  for  this  program. 

Nominally,  the  target  VLMA  fiber  has  a  SOpm  effective  core  diameter  with 
stable,  diffraction- limited  emission.  We  have  developed  several  guidelines  motivating 
our  concepts  and  designs,  which  will  be  discussed  below: 

(a)  Significant  HOM  filtering  -  If  HOMs  have  high  loss  relative  to  the  fundamental 
mode,  the  instability  threshold  will  be  higher  and  stable  beam  quality  will  result. 
In  other  words,  power  seeded  from  the  LPoi  mode  to  an  HOM  that  experiences 
loss  results  in  lower  net  effective  mode  coupling.  The  utility  of  this  idea  with 
respect  to  modal  instabilities  was  experimentally  proven  by  ARFL,  and  later 
numerically  corroborated  by  AS  Photonics.  The  first  guideline  is  therefore  to 
design  waveguides  with  strong  mode  filtering  of  some  kind.  This  is  also  viewed 
as  the  most  challenging,  since  to  date  no  such  waveguides  exist  commercially 
beyond  25  pm  core  diameter. 

(b)  Large  mode-index  contrast  -While  the  LPn  mode  is  typically  the  most 
detrimental  to  LPoi  propagation  in  general,  the  small  modal  index  difference 
between  these  modes  leads  to  a  large  beat  length  and  correspondingly  significant 
phase  matching,  which  are  ripe  conditions  for  initiating  modal  instabilities.  By 
designing  for  larger  spacing  of  effective  indices  between  the  LPoi  and  the  HOMs, 
the  beat  length  will  be  shorter,  reducing  the  phase  matching.  Alternatively 
viewed  from  the  grating  perspective,  longer  period  gratings  necessarily  have 
fewer  “reflections”  per  unit  length  that  add  coherently,  necessitating  a  longer  or 
stronger  grating  to  generate  the  same  effect.  In  addition,  the  longitudinal  heat 
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flow  (which  we  have  experimentally  measured  as  significant  over  cm-seale 
lengths  in  some  eases  [24])  should  aid  in  diffusing  shorter  period  gratings.  The 
seeond  guide  is  therefore  to  maximize  the  differential  modal  index  between  the 
fundamental  mode  and  the  pertinent  HOMs  (which  may  not  be  the  LPn 
depending  on  the  speeific  fiber  design). 

(e)  Redueed  spatial  overlap  -The  perturbative  eoupling  that  leads  to  modal 
instability  depends  on  the  spatial  intensity  overlap  between  the  LPoi  and  HOMs  in 
the  region  of  the  (index)  perturbation.  By  designing  waveguides  that  suffieiently 
perturb  the  spatial  profiles  from  the  nominal  (step-index)  LP  mode  profiles,  the 
spatial  overlap  and  thus  coupling  between  modes  ean  be  reduced,  increasing  the 
instability  threshold.  In  addition,  this  joint  mode  overlap  is  only  relevant  in  the 
regions  of  thermal  perturbations,  which  are  dictated  by  the  gain  distribution.  The 
third  guide  is  therefore  to  reduce  the  spatial  (intensity)  overlap  between  the  LPoi 
and  the  relevant  HOMs  partieularly  over  the  region  that  is  doped  with  rare-earth 
ions  (e.g.,  ytterbium  in  the  ease  of  interest  to  this  program). 

(d)  Redueed  mode  set  -  Nominally,  the  modal  effeetive  indiees  (whose  differenee 
govern  the  mode  eoupling)  are  bound  between  the  core  and  cladding  refractive 
indiees.  If  designs  can  be  reached  that  reduce  the  total  number  of  modes  within 
the  same  eigenspaee  (ncore  >  ngff  >  ndad),  then  it  will  likely  also  inerease  the 
differential  modal  index  [as  desired  per  guideline  (b)].  In  addition,  it  will  also 
reduee  the  number  of  degrees  of  freedom  in  the  system,  whieh  should 
fundamentally  lead  to  inereased  stability  sinee  there  is  less  “noise”  to  seed  the 
instability  process  (e.g.,  from  imperfeet  seed  launeh).  The  fourth  guide  is 
therefore  to  reduce  the  total  number  of  modes  within  a  given  NA  space. 

Note  that  nature  provides  various  ways  for  several  of  these  guidelines  to  work 
together.  For  example,  modes  that  do  not  have  significant  spatial  overlap  [guideline  (c)] 
are  likely  to  have  larger  effeetive  index  contrast  [guideline  (b)]  sinee  the  effective  index 
is  related  to  the  spatial  overlap  between  the  optieal  field  and  the  refraetive  index.  In 
another  example,  using  strong  gain  filtering  (50-60%  radial  gain  eonfmement)  to  eut  out 
the  LPii  mode  [guideline  (a)]  means  that  the  next  HOM  will  eause  the  problems; 
however,  the  next  HOM  will  have  a  lower  effeetive  index  [guideline  (b)]  and  redueed 
spatial  overlap  [guideline  (c)]. 

Many  sueh  eandidates  were  eoneeptualized,  modeled  for  spatial  overlap  and  Angff, 
and  assessed  using  the  REM  model.  The  key  suceessful  design  results  are  described  in 
the  following  subseetions. 

6.1  Confined-Gain  Fiber 

Confined  gain  has  been  recognized  as  one  of  the  key  methods  for  mode  eontrol  in  VLMA 
fibers,  therefore  conforming  to  guideline  (a).  Not  only  has  it  been  well  modeled  [6,8], 
but  it  has  been  experimentally  validated  [25]  ineluding  powers  up  to  300  W  [26].  The 
key  to  “gain  filtering”  is  not  to  introduce  loss  into  the  amplifier  but  rather  to  deny  gain  to 
HOMs  under  all  eircumstances.  The  physics  of  gain  filtering  using  confined  gain  is  best 
illustrated  in  Fig.  11. 

Although  the  impaet  of  eonfmed  gain  on  the  TMI  threshold  was  originally  studied 
in  the  AFRL  publications  [15],  the  radial  confinement  was  limited  to  80%  of  the  eore 
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diameter.  While  the  results  still  showed  an  inerease  in  the  TMI,  the  full  potential  of 
eonfmed  gain  was  not  realized  in  those  simulations.  Theoretieal  studies  [6,8]  report  that 
50-60%  of  the  eore  diameter  represents  the  best  modal  suppression.  Sinee  filtering  out 
ROMs  neeessarily  results  in  a  reduetion  of  TMI  tendeneies  (via  inereased  threshold),  it  is 
prudent  to  optimize  the  effeetiveness  of  gain  filtering. 


Fundamental 
mode  launched 


Gain  is  spatially 
saturated 


HOMs  overlap  with 
remaining  gain 


Confined  gain 
(centrally  doped) 


Transverse  coordinate 


Figure  11:  Conventional  gain  fibers  (top)  provide  gain  to  HOMs  especially  when  the  gain  is 
saturated.  Confined  gain  (bottom)  eliminates  gain  for  HOMs  even  under  saturation  conditions. 


In  our  simulations,  we  used  a  50|j,m  core  diameter  with  a  60%  doping  diameter. 
With  the  same  cladding  diameter,  this  represents  the  same  approximate  laser  kinetics  as  a 
SOpm  LMA  fiber  core.  It  is  important  to  note  that  confining  the  gain  does  not  mean 
increasing  the  absorption  length  if  the  waveguide  diameter  is  increased. 

The  results  of  the  simulations  were  extremely  encouraging.  While  stability  was 
maintained  up  to  IkW  of  power,  higher  powers  also  resulted  in  another  stability  region. 
This  means  that  at  the  very  minimum,  over  a  factor  of  two  was  gained  in  the  TMI 
threshold  for  a  SOpm  fiber.  It  is  notable  that  the  TMI  threshold  for  the  50pm  confined- 
gain  fiber  is  at  least  40%  higher  than  the  TMI  threshold  for  a  conventional  25pm  LMA 
fiber.  These  results  are  captured  in  Fig.  12,  where  the  confined-gain  simulation  shows 
both  observed  thresholds. 

These  extremely  encouraging  results  necessitate  further  study  of  confined-gain 
fiber  as  a  new  platform  for  VLMA  fiber  types.  In  addition  to  continued  power  scaling, 
more  modes  need  to  be  included  in  the  study;  since  the  LPn  mode  is  typically  eliminated 
by  gain  filtering  due  to  the  null  in  the  center  of  its  intensity  profile,  the  LP02  mode 
becomes  the  next  threat  to  beam  quality  and  mode  stability. 
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Figure  12:  TMI  threshold  for  confined-gain  fiber  (orange)  compared  to  conventional  LMA  fiber 
(blue).  The  lines  are  guides  for  the  eye,  and  the  error  bars  represent  the  pump  increment.  The 
dual  points  for  the  confined-gain  fiber  are  explained  in  the  text. 

6.2  Cladded  Linear  Index  Graded  (CLING)  Fiber 

In  addition  to  gain  tailoring,  the  concept  of  refractive-index  profiling  was  also  brought  to 
bear  in  this  program.  In  the  simplest  terms,  waveguide  modes  are  well  represented  by  a 
quantum  mechanical  potential  (i.e.,  a  quantum  well)  viewed  upside  down.  The 
eigenvalues  (effective  indices)  are  positioned  in  the  “well”  at  an  appropriate  “energy 
level”  (location  in  refractive  index  space),  and  the  mode  size  is  largely  dictated  by  the 
location  of  the  edges  of  the  “well”  (refractive  index  boundary).  This  concept  is  shown  in 
Fig.  13(a).  Note  that  the  modes  in  this  case  are  approximately  the  same  size  and  thus 
have  high  spatial  (intensity)  overlap. 

One  method  to  reduce  the  overlap  between  the  modes,  in  order  to  follow 
guideline  (c),  is  to  alter  the  shape  of  the  “well”.  Based  on  the  above  logic,  it  is 
reasonable  to  assume  that  if  the  “well”  were  tapered  (i.e.,  a  graded-index),  the  second 
mode  would  be  larger  than  the  fundamental  mode.  This  concept  is  depicted  in  Fig.  13(b), 
and  is  oversimplified  since  the  changes  in  ngff  were  not  taken  into  account.  Nonetheless, 
the  concept  of  modifying  the  waveguide  profile  to  affect  the  mode  sizes  stands. 

In  studying  this  concept,  the  shape  of  the  index  profile  was  found  to  be  important. 
Although  only  two  modes  are  shown  in  Fig.  13,  a  VLMA  fiber  will  by  default  have  many 
modes,  with  the  effective  indices  of  the  lowest  modes  gathered  higher  in  the  inverted 
index  well.  Designing  a  conventional  parabolic  graded  index  provided  challenges  since 
the  relative  LPoi-LPn  overlap  was  strongly  dependent  the  precision  to  which  the 
parameters  were  hit;  a  5%  change  in  the  width  of  the  waveguide  resulted  in  a  significant 
reduction  in  the  desired  size  mismatch.  With  a  linear  index  grade,  the  sensitivity  to 
manufacturing  tolerances  was  found  to  be  significantly  reduced;  if  both  eigenvalues 
moved  up  or  down  the  well,  they  largely  did  so  together  with  a  reasonably  consistent 
overlap.  A  sub-linear  index  profile  did  not  provide  sufficient  mode  separation,  much  like 
the  steep  edges  of  a  conventional  step-index  fiber. 
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Figure  13:  Comparison  of  modes  between  step-index  and  graded-index  fibers,  (a)  Modes  of  a 
conventional  step-index  fiber  are  nearly  the  same  size,  (b)  In  contrast,  modes  of  a  graded-index 
fiber  are  not  the  same  size,  with  higher-order  modes  being  increasingly  larger. 

Therefore,  due  to  both  physies  and  toleranees,  we  settled  on  the  cladded  linear 
index  graded  fiber,  or  CLING  fiber.  The  speeific  design  we  selected  is  shown  in  Fig.  14. 
The  initial  idea  was  to  modify  a  conventional  step-index  fiber  with  a  0.06  NA  by 
bumping  the  center  up  and  the  edges  down  by  the  same  amount.  The  total  core  diameter 
was  expanded  to  100  pm  in  order  to  allow  sufficient  expansion  of  the  LPn  mode. 


Figure  14:  Modes  of  a  LING  fiber  compared  to  modes  of  a  conventional  50-pm  step-index  fiber 
of  the  same  average  NA  (0.06).  The  edge  index  step  of  the  LING  fiber  is  0.0041  while  the  graded 
portion  spans  0.0011. 

Figure  14  shows  both  the  LPoi  and  LPn  modes  for  two  fibers  having  the  same 
effective  area  for  the  fundamental  mode:  the  CLING  fiber  (blue  and  orange,  respectively) 
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and  a  50)j,m  LMA  fiber  (purple  and  green,  respeetively).  Note  that  although  the 
fundamental  modes  are  approximately  the  same  size  and  shape,  the  LPn  mode  of  the 
LING  fiber  is  substantially  (almost  40%)  wider.  Radially,  this  means  that  the  2D  overlap 
is  redueed  by  nearly  a  factor  of  two,  in  principle  significantly  raising  the  TMI  threshold. 

Figure  15  shows  the  TMI  threshold  of  the  CLING  fiber  in  comparison  to  the 
conventional  step-index  LMA  fiber  type.  The  CLING  fiber  shows  a  68%  increase  in  the 
TMI  laser  threshold  over  a  conventional  LMA  fiber  of  the  same  fundamental  mode  size 
(50pm  effective  core  diameter).  Moreover,  the  TMI  threshold  is  higher  than  that  of  the 
industry  standard  25 jam  LMA  fiber,  but  with  a  mode  area  that  is  4x  larger. 


20  25  30  35  40  45  50 

Fiber  Core  Diameter  (pm) 


Figure  15:  TMI  threshold  for  CLING  fiber  (purple)  compared  to  conventional  (step-index)  LMA 
fiber  (blue).  The  lines  are  guides  for  the  eye,  and  the  error  bars  represent  the  pump  increment. 

The  CLING  fiber  modeled  here  had  a  50pm  gain-dopant  diameter  in  order  to  be 
similar  to  conventional  LMA  fiber.  In  spite  of  the  difference  in  LPqi  and  LPn  mode 
sizes,  the  concept  of  gain  filtering  has  not  been  fiilly  applied  here.  Adding  gain  filtering 
to  this  concept  would  require  a  gain  doping  diameter  of  ~30pm,  much  smaller  than  what 
was  used  in  the  present  simulations.  Additional  exploration  and  design  of  CLING  fibers 
is  highly  recommended,  particularly  with  the  inclusion  of  gain  filtering. 

6.3  Trefoil  Fiber 

In  an  effort  to  meet  both  guideline  (c)  (reduced  overlap)  and  guideline  (d)  (reduced  mode 
set),  a  triangular  core  was  considered.  Although  such  a  core  could  in  principle  be 
fabricated,  exploiting  conventional  rod-in-tube  or  stacked-rod  geometries  are  much  more 
practical.  Therefore,  the  geometry  considered  is  that  of  three  rods  in  what  we  call  the 
Trefoil  fiber.  Each  of  the  three  cores  has  0.06  NA,  with  a  diameter  selected  to  produce  an 
effective  area  of  the  fundamental  mode  that  is  equivalent  to  a  50  pm  step-index  fiber  of 
the  same  NA. 

The  first  three  modes  of  the  Trefoil  fiber  investigated  are  shown  in  Fig.  16.  In 
this  figure,  the  green  lines  represent  the  core  boundaries,  while  the  red  ovals  denote  the 
location  of  power  in  the  fundamental  mode.  One  immediate  success  of  this  fiber  type  is 
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that  the  LPn^a  mode  has  at  most  a  60%  spatial  overlap  with  the  LPoi  mode.  Perhaps  more 
importantly,  the  degeneracy  of  the  LPn  modes  has  been  broken.  The  LPn^b  mode  has  a 
much  lower  overlap  with  the  fundamental  mode,  near  30%.  This  means  that  in  practice, 
only  half  of  the  power  will  be  leaked  out  of  the  fundamental  mode,  which  should  result  in 
not  only  better  beam  quality,  but  higher  TMI  threshold. 


Figure  16:  Lowest-order  modes  of  a  Trefoil  fiber.  The  green  circles  annotate  the  core  locations, 
while  the  red  ovals  are  guides  to  the  eye  for  marking  the  fundamental  mode  power. 

Unfortunately,  since  the  concept  came  late  in  the  program,  only  a  single  set  of 
TMI  simulations  were  performed  for  this  preliminary  Trefoil  fiber  design,  with  the  results 
shown  in  Fig.  17.  The  TMI  threshold  of  650W  showed  a  promising  37%  increase  over 
the  50pm  LMA  fiber  even  for  this  un-optimized  design.  However,  regions  of  stability 
were  found  at  pump  powers  up  to  875W.  Given  the  reduced  spatial  overlap  and  the 
multiple  regions  of  stability,  additional  exploration  and  design  of  Trefoil  fibers  is  highly 
recommended,  particularly  with  the  inclusion  of  gain  filtering  and  attention  paid  to  the 
beam  quality  of  the  fundamental  mode. 
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Figure  17:  TMI  threshold  for  Trefoil  fiber  (teal)  compared  to  conventional  (step-index)  LMA 
fiber  (blue).  The  lines  are  guides  for  the  eye,  and  the  error  bars  represent  the  pump  increment. 
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7.  Summary  and  Conclusions 

Many  different  fiber  types  and  designs  were  explored  in  this  program  using  a  validated 
rate-equation-based  model  for  accurately  modeling  the  beam  quality  in  VLMA  fiber 
amplifier  while  simultaneously  including  the  effects  of  thermal  lensing  and  thermal  mode 
instability. 

Mode-area  scaling  using  conventional  LMA  fibers  is  simply  not  possible  since  the 
TMI  threshold  decreases  with  increased  core  area.  Modifications  in  pump  geometry  can 
provide  a  modest  benefit,  particularly  when  flattening  out  the  temperature  profile  along 
the  fiber 

Confined  gain  was  studied  and  found  to  provide  exceptional  resilience  to  TMI  (at 
least  2x  increase  in  TMI  threshold),  as  was  originally  proposed.  Two  new  promising 
fiber  types,  CLING  and  Trefoil,  were  invented,  and  both  exhibit  features  with  modeled 
and  expected  resilience  to  TMI. 


Figure  18:  TMI  threshold  results  summary  for  all  conditions  and  fiber  types.  TMI  threshold  is 
shown  for  forward  (blue),  backward  (red),  and  bi-directional  (green)  pumping  as  a  function  of 
core  diameter.  Also  shown  are  the  TMI  threshold  for  the  confined-gain  fiber  (orange),  the 
CLING  fiber  (purple),  and  the  Trefoil  fiber  (teal).  Note  that  the  regions  of  stability  at  higher 
pump  powers  were  found  for  the  confined-gain  and  Trefoil  fiber  types.  The  lines  are  guides  for 
the  eye,  and  the  error  bars  represent  the  pump  increment. 

The  most  important  results  are  captured  in  Fig.  18,  which  shows  all  of  the 
calculated  TMI  thresholds  together  in  a  single  plot,  and  Table  1,  which  quantifies  and 
annotates  the  findings  specifically  for  scaling  to  50  pm  core  diameter.  It  should  also  be 
noted  that  multiple  effects  can  be  integrated  into  a  single  design,  such  as  using  the 
CLING  fiber  with  gain  filtering  (confined  gain)  and  bi-directional  pumping,  for 
maximum  TMI  suppression. 
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Table  I:  Program  summary  results  for  50/am  effective  core  diameter  fibers 


Scaling  method 

TMI  Threshold 

Notes 

LMA  core  scaling 

37%  decrease 

Not  viable  -  new  fiber  needed 

Bi-directional  pumping 


No  change 


10%  increase 


>110%  increase 


68%  increase 


37%  increase 


Optimize  to  minimize  temperature  peak 
Even  higher  power  stability  regions 
No  explicit  gain  filtering  included 
Un-optimized  design,  no  gain  filtering 


Key  results  from  this  program  will  be  presented  at  conferences  and  written  up  for 
publication  in  refereed  journals.  The  program  results  will  be  summarized  in  a 
presentation  at  the  HEL-JTO  review  during  the  week  of  May. 

8.  List  of  Symbols,  Abbreviations,  and  Acronyms 


ADI 

alternating  direction  implicit 

AFRL 

Air  Force  Research  Laboratory 

BPM 

beam  propagation  method 

CLING 

cladded  linear  index  graded 

CW 

continuous  wave 

DEW 

directed  energy  weapon 

FD-BPM 

finite  difference  beam  propagation  method 

FEA 

finite  element  analysis 

HOM 

higher  order  mode 

LMA 

large  mode  area 

NA 

numerical  aperture 

REM 

rate  equation  modeling 

TMI 

thermal  mode  instability 

VLMA 

very  large  mode  area 
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